The doubly heavies: 
QQqq and QQqq tetraquarks and QQq baryons 
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>• . Abstract 

IT) ; 

IT ■ Recent discoveries by Belle and BESIII of charged exotic quarkonium- 

^^ ! like resonances provide fresh impetus for study of heavy exotic hadrons. 

^ | In the limit N c — > oo, Mq — > oo the QQqq' tetraquarks (TQ-s) are 

O ■ expected to be narrow and slightly below or above the Qq' and Qq 

two-meson threshold. The isoscalar TQ-s manifest themselves by de- 

• • | cay to (QQ) rnr, and the ~30 MeV heavier charged isotriplet TQ-s 

. -h ■ by decays into (QQ)ir. The new data strongly suggest that the real 

^% . world with N c = 3, Q = c,b and q, q' = u,d is qualitatively described 

j^ I by the above limit. We discuss the relevant theoretical estimates and 

suggest new signatures for TQ-s in light of the recent discoveries. We 
also consider "baryon-like" states QQ'qq', which if found will be di- 
rect evidence not just for near-threshold binding of two heavy mesons, 
but for genuine tetraquarks with novel color networks. We stress the 
importance of experimental search for doubly-heavy baryons in this 
context. 
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It has been realized early on that quark models and QCD sustain a much 
richer pattern of different multi-quark and/or color network configurations, 
beyond the "non-exotic" standard qq mesons and qqq baryons. Still, produc- 
tion rates of such particles are often suppressed and the light pions will in 
most cases allow rapid decays of the exotics into final states with pion(s) - 
turning them into very broad resonances. This explains why the vast ma- 
jority of known hadrons [TJ are just simple mesons and baryons. There is, 
however, growing evidence for the existence of exotic mesons containing both 
heavy and light quark- ant iquark pairs - ccqq [2] and as well as manifestly 
exotic bbdu [3], [I] and most recently ccdu [5]. 

In a recent article [6] S. Weinberg made prescient observations on qqqq 
tetraquarks (TQ-s) in the large- N c limit, correcting some statements in S. Cole- 
man's Erice Lectures j7]. Specifically, TQ-s can be created from the vacuum, 
as in lattice calculations [Sj, in two ways, via products of color-singlet qq 
bilinears: 

|T) = [gi(aOlWx)] ■ [q 3 (x)Tjq 4 (x)]\0) (1) 

and the "crossed" 2 <^> 4 version, with T indicating a Dirac algebra matrix. 
TQ-s with a specific value of isospin can be constructed by taking linear com- 
binations of the {1, 2} or {3, 4} terms with their {2, 1} or {4, 3} permutation, 
respectively. 

The rapid two-pion decay of /o(980) and rj7i decay of ao(980) are avoided 
if, (as first suggested in a bag model context by Jaffe [9]), the latter are 
ss(uu ± dd) TQ-s. The above decays are then slowed by Zweig's rule, 1/N 2 
suppression of ss — > uu/dd transitions, and the "natural" decay channel KK 
is above 980 MeV. 

These features stem from the fact that the strange quark is somewhat 
heavier than the light u and d quarks, making the / (980) and a (980) the 
harbingers of recent and likely future TQ-s, QQqq with Q = c,b; q = u,d. 
Such systems have an additional relevant small parameter, namely m q /MQ ~ 
1/5 and ~ 1/15, respectively. Rearrangement of QQqq subsystems with sub- 
sequent decay into J/ip (or T) and pion(s) is always allowed by kinematics, 
even when the TQ is lighter than the sum of the masses of the Qq + Qq 
mesons, but such rearrangements are 1/N C suppressed. 

This suppression is the key to the relatively small width of the manifestly 
exotic charged resonances recently discovered in both bb by Belle [I] - and 
just a few days ago, the remarkable peak Z c (3900) at 3899.0 ± 3.6 ± 4.9 MeV 



with T = 46 ± 10 ± 20 MeV in cc system by BESIII [5J J3 

The relatively slow decay of these exotic resonances implies that the dom- 
inant configuration of the QQqq four-body system is not that of a low-lying 
QQ quarkonium and pion(s). The latter have a much lower energy than 
the respective two-meson thresholds MM* and M*M*, (M = D,B), but do 
readily fall apart into (QQ) and pion(s) and would result in very large decay 
widths. 

Instead, just like in the case of / (980) and a (980), we should view these 
systems as loosely bound states and/or near threshold resonances in the two 
heavy-meson system. 

Such "molecular" states, DD*, etc., were introduced [TTJ Q2] in analogy 
with the deuteron which binds via exchange of pions and other light mesons, 
and were referred to as "demons" . The key observation is that the coupling 
to the heavy mesons of the light mesons exchanged (n, p, etc.) becomes 
universal and independent of Mq for Mq — > oo, and so does the resulting 
potential in any given J CP and isospin channel. In this limit the kinetic 
energy ~ p 2 /(Mq) vanishes, and the two heavy mesons bind with a binding 
energy ~ the maximal depth of the attractive meson-exchange potential. 

For a long time it was an important open question whether these con- 
sideration apply in the real world with large but finite masses of the D and 
B mesons. The recent experimental results of Belle [I] and BESIII |5J, to- 
gether with theoretical analysis in Refs. [19] and [20] strongly indicate that 
such exotic states do exist - some were already found and more are predicted 
below. 

Due to parity conservation, the pion cannot be exchanged in the MM 
system, but it does contribute in the MM* and M*M* channels. The t\ ■ t s 
isospin nature of the exchange implies that the the binding is 3 times stronger 
in the isoscalar channel, experimentally accessed in the decay into Tttti or 
J/i/j7in, than in the isovector channel, seen in decay into T7r or J/ipir. It 
was estimated JT9J |20] that in the bottomonium system this difference in 
the binding potentials raises the 1=1 exotics by up to 40-50 MeV above 



'''Note that a O(30) MeV energy implies a velocity of ^(1/8) c for D and D* in 
the final state. For a system of size O(0.7) Fermi, this implies a traversal time and the 
corresponding crude estimate for the S'-wave state width of 35 MeV. The fact that this is 
close to the observed width suggests that we indeed have a loosely-bound DD* resonance 
decaying mainly into its constituent mesons. 



the 1=0 exotics. In the charmonium system this splitting is expected to be 
slightly smaller, because the DD*/D*D* states are larger than BB* / B*B* 
and therefore the potential in both 1=0 and 1=1 channels is expected to be 

somewhat weaker. 

We are interested in the dependence of the binding energy on the heavy 
meson masses in the MM, MM* and M*M* systems controlled by Hamil- 
tonians of the form: 

H = a-p 2 + V{r) (2) 

where a = l/2fi Icd , with fi TeA the reduced mass and the channel-specific po- 
tential V(r) is assumed to be independent of the heavy-meson mass. 

For Hamiltonians linear in a parameter a, the ground-state energy Eq(o) 
is a convex function of a, since second-order perturbation makes a nega- 
tive contribution to the ground state. For H = J27 a iHi, the function 
Eo(ai, a 2 , . . . , a n ) describes a convex surface in n dimensions. With addi- 
tional theoretical inputs or experimental results, the n-dimensional convexity 
might help estimating how the binding is affected by scaling various pieces 
of the potential. 

Convexity implies that Eq(o) at any point in the (0, A) interval exceeds 
the linear interpolation 

E (a) > E (0) + a (3) 

Related arguments were used to prove flavor-dependent baryon-meson mass 
inequalities [21, 22J. Convexity was also used to predict the masses of doubly- 
heavy be mesons [23] . 

Convexity is not affected by linear transformations of the a, and we can 
use any unit of inverse mass to convert the parameter a ~ 1/Mq into a 
dimensionless parameter a. 

Since we are comparing with the pion and with other light exchanges, we 
use the "natural" unit of ~ 0.7 Fermi ~ 3.5 GeV -1 for the range and width 
of the corresponding potentials. 

With mr> ~ 2 GeV and tub ~ 5.3 GeV the relevant dimensionless pa- 
rameters a(D) = 1/7 and a(B) = 1/18.5 are indeed small, suggesting that 
even the first-order perturbation linear interpolation may be satisfactory. 



In the infinite mass case, corresponding to a = 0, the relative distance 
between the mesons is effectively fixed at the minimum of the potential and 
the binding Eq(0) is then just the attractive V m tn at this position. 

We can use the existing data in order to make a very rough estimate of 
the isovector binding potential which is E\ =x in the ttiq — > oo limit, i.e. at 
a = 0. We have two data points: Z c (3900) at a(D) is approximately 27 MeV 
above DD* threshold and ^(10610) at a(B) is approximately 3 MeV above 
BB* threshold. Linear extrapolation to a = yields El =l (a=0) fa —11.6 
MeV. In view of the convexity, the actual binding energy is likely to slightly 
exceed this linear extrapolation. 

We can then use this result for the isovector channel to estimate the 
BB* binding in the isoscalar channel. Assuming that the isoscalar binding 
energy in the iuq —¥ oo limit is 3 times larger than for the isovector, i.e. 
El= Q (a=ti) fa 3 ■ (-11.6) = -35 MeV. V(3872) is at DD* threshold, pro- 
viding an additional data point of E[ =0 (a(D)) m in the isoscalar channel^ 
Linear extrapolation to a(B) yields approximately —20 MeV as the BB* 
binding energy in the isoscalar channel. 

The upshot is that the newly discovered Z c (3900) isovector resonance 
confirms and refines the estimates in [191 120] for the mass of the putative 
BB* isoscalar bound state. This immediately leads to several predictions: 

a) two 1 = narrow resonances in the bottomonium system, about 23 
MeV below Z 6 (106010) and Z 6 (10650), i.e. about 20 MeV below the 
corresponding BB* and B*B* thresholds; 

b) an / = 1 resonance about 30 MeV above D*D* threshold; 

c) an / = resonance near D*D* threshold. 



Various authors, e.g. [TT],jT2] considered QQq\q2 TQ-s. If such states do 
exist, producing and discovering even the lightest ccud is an extraordinary 
challenge. One needs to produce two QQ pairs and then rearrange them, so as 
to form QQ and QQ diquarks, rather than the more favorable configuration 
of two QQ and color singlets. Then the QQ diquark needs to pick up a 



■^For simplicity we assume X(3872) is an isoscalar, since it has no charged partners, 
and we ignore here the issue of isospin breaking in its decays. 



ud light diquark, rather than a q, to make a QQq baryon, suppressing the 
production rate of these TQ-s below the rate of QQq production. 

A small ray of hope comes from the observation of the doubly-heavy 
B c = (be) mesons [13], suggesting that simultaneous production of bb and cc 
pairs which are close to each other in space and in rapidity and can coalesce 
to form doubly-heavy hadrons is not too rare. For example, in the last paper 
in Ref. JT3J CDF reported 108 ± 15 candidate events in the Bf — >• J/ijnr* 1 
channel. This is an encouraging sign for the prospects of producing and 
observing the ccq and bes baryons and hopefully also the ccud TQ. ATLAS 
and CMS and especially LHCb probably have the best chance of discovering 
these states. 

If the new TQ lies below say, the DD* threshold, it will be stable under 
the strong interaction and will decay only weaklyj^l As we discuss below, if 
they lay above threshold, they may still may manifest themselves as narrow 
DD* resonances. 

The following theoretical considerations are relevant for the last point 
above. First, in clear contradistinction to the states discussed earlier, the 
new TQ-s are unlikely to be molecular. This is because both heavy mesons 
contain light antiquarks, rather than a q in one and a q in the other, causing 
the uj and p exchange to be repulsive, rather than attractive. With no bound 
DD* states, one might expect the new DD* states to lay higher and manifest 
themselves as broad resonances, at best. 

This bleak picture will change dramatically if the QQud TQ is not a 
simple molecular state, but instead a connected color network, consisting of 
a QQ diquark coupled to a 3*, with a matching ud diquark coupled to a 
3 C . The QQ diquark will then be bound by the large Coulombic interaction 
0(o%Mq), as discussed below. Finding such a novel color network would be 
particularly exciting [12], [16]. With some poetic license, we can compare it 
to the discovery of the Cqq buckyballs. 

The Coulombic binding of the QQ diquark, 

E b (QQ) = i ■ 1 ■ \a\M Q = l -a 2 s M Q (4) 



§A statement slightly modified if the TQ mass is above the 2mc threshold, as it can 
then first decay through an exceedingly narrow EM decay D* — > Z?°7 [Tj,[T4], followed by 
the weak decays of the the D-s [12] • 



is half as strong as that of a color-singlet QQ. Another | factor reflects 
the need to use the reduced mass, and the \ol\Mq is the standard hydro- 
gen binding formula. Thus the heavy diquark QQ becomes infinitely deeply 
bound as Mq — y oo. For m c ~ 1.6 GeV and using 1/2 for the relevant a s , 
the Coulombic binding ~ 50 MeV is moderate. However, the | ratio of QCD 
Coulomb interaction in the diquark, versus its strength in the QQ color sin- 
glet system likely extends to all QCD interactions, including the confining 
stringy potential [13] . 

If indeed the ccud TQ is mostly described by the new color network, then 
its decay into DD* will be suppressed by the required rearrangement of the 
color. 

In view of the above discussion, we now focus on the question whether 
the lowest state of the new color network is below or near the DD* threshold. 

There are interesting parallels between the QQqq TQ-s and doubly-heavy 
baryons QQq. In both types of systems there is a light color triplet - a quark 
or an anti-diquark - bound to a heavy diquark. Because of this similarity, 
experimental observation of doubly-heavy baryons is very important not just 
in its own right, but as source of extremely valuable information for deducing 
the properties of the more exotic QQqq tetraquarks. Such deduction can 
carried out just as it was done for 6-baryons. 

In the last few years it became possible to accurately predict at the level 
of 2-3 MeV the masses of heavy baryons containing the 6-quark: T,b(bqq), 
Eb(bsq) and fib(bss) [2H (23J [26] . These predictions used as input the masses 
of the B, B s , D and D s mesons, together with the masses of the corresponding 
c-baryons T, c (cqq), E c (csq) and Q c (css). 

An analogous approach to the masses of the QQqq, suggested the relation 
m(ccud) = m(E ccu ) + m(A c ) - m(D°) - ~[m(D*) - m(D)] (5) 



designed to optimally match interactions on both sides [12J. To date, only 
the SELEX experiment at Fermilab reported doubly charmed ccd and ecu 
baryons with mass ~ 3520 MeV [27] - a result not confirmed by other ex- 
periments [281 [29] . Substituting this mass in eq. (J3J) as a placeholder in a 
"proof of concept" estimate yields M(ccud) ~ 3900 MeV, which is just 30 
MeV above the D*D threshold. One should keep in mind that, as implied by 
eq. (jlj), the binding in (bbqq) is expected to be significantly stronger than in 



(ccqq). In any case, reliable experimental information on QQq baryon masses 
is clearly essential for settling the issue. 



Note added 

While this work was being finalized, we became aware of several papers deal- 
ing with closely-related subjects, using different approaches 
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